Connection between charge transfer and alloying core-level shifts based on 

density-functional calculations 



M. Methfessel 

Institute for Semiconductor Physics (IHP), W alter- K or sing- Str. 2, D-15230 Frankfurt (Oder), Germany 

Vincenzo Fiorentini and Sabrina Oppo 

Istituto Nazionale per la Fisica della Materia - Dipartimento di Fisica, Universita di Cagliari, Italy 
(to appear on Phys. Rev. B, 15-February-2000) 

The measurement of alloying core-level binding energy (CLBE) shifts has been used to give a 
precise meaning to the fundamental concept of charge transfer. Here, ab-initio density-functional 
calculations for the intermetallic compound MgAu are used to investigate models which try to 
make a connection between the core levels shifts and charge transfer. The calculated CLBE shifts 
agree well with experiment, and permit an unambiguous separation into initial-state and screening 
contributions. Interestingly, the screening contribution is large and cannot be neglected in any 
reasonable description. Comparison of the calculated results with the predictions of simple models 
show that these models are not adequate to describe the realistic situation. On the positive side, 
the accuracy of the density-functional calculations indicates that the combination of experiments 
with such calculations is a powerful tool to investigate unknown systems. 
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I. INTRODUCTION 

The concept of charge transfer is fundamental to chem- 
istry and condensed-matter physics. Unfortunately, it is 
frustratingly difficult to give a precise definition of charge 
transfer, or even a well-defined prescription for measur- 
ing it. This is equally true for related quantities such as 
electronegativity and bond ionicity. 

There have been a number of attempts to relate the 
charge transfer in an alloy to the positions of the core 
levels. These energies can be measured with high ac- 
curacy using X-ray spectroscopy!!! In general terms, the 
core levels of an atom are shifted when the atom's en- 
vironment changes. Interesting cases are, fox example, 
when a crystal is formed, out of free atoms ,□ when an 
atom is at the surfaceETtl rather than in the bulk of a 
solid, gp. when an alloy is formed out of two elemental 
solids .LrU 

In the case of the alloy core-level shift (the subject of 
this paper) a major objective has been to find a well- 
defined connection between the measured shift and the 
charge transfer between the constituents. Clearly, when 
charge is moved from one atom to another, an electro- 
static potential builds up, which modifies the energy 
needed to eject an electron from a core level into the vac- 
uum. In the simplest form of the "potential model" , the 
change of the potential felt by a core electron is described 
using a M adelung term and an on-site contribution.!^ Un- 
fortunately, cancellation between these effects and uncer- 
tainty in the the model parameters make it difficult to 
extract reliable charge transfers using this approach. 

Furthermore, the simple potential model is valid only 
for the "initial state" picture, i.e. when describing the 
positions of the core levels in the alloy and the pure metal 



before a core electron is removed. To compare with the 
measured binding energies, a final-state screening contri- 
bution must be taken into account: after a core hole is 
created, the remaining electrons relax to screen the hole. 
The kinetic energy of the emitted electron includes the 
screening energy. This can be included in the formalism 
(by a term generally denoted 6R), but this adds yet an- 
other parameter whose numerical value is poorly known. 

An alternative procedure is described in Ref. |8[ It 
was pointed out that the relaxation energy can be mea- 
sured directly via the shift of the Auger parameter (the 
sum of core- level ionization and Auger-energy shifts). For 
the compounds AuMg and AuZn, the resulting values of 
5R axe inconsistent with the estimates used in an earlier 
worMO based on the potential model, even though plau- 
sible ionicities were deduced there. A modified potential 
model was presentedjj which relates the valence charge 
transfer in a metal to an atomic property, namely the 
change in the potential at the core due to changes in the 
valence and core occupation numbers. 

Under these circumstances, it is desirable to make a 
detailed theoretical analysis of a typical system using a 
method which can quantify the various contributions un- 
ambiguously. Herep we use ab-initio density-functional 
total-energy calculations to study the MgAu alloy. Us- 
ing a supercell technique, the Mg, Au, and MgAu metals 
with and without a core hole on a selected atom can 
be described accurately. This cleanly separates the core- 
level shift into initial-state and final-state relaxation con- 
tributions, which can then be checked against the appro- 
priate models. In addition, direct inspection of the densi- 
ties of states, on-site charges, and screening charge distri- 
butions gives an understanding of the effects of alloying 
and the different screening responses to a core hole. 

/From the results, we are led to the conclusion that 
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the final-state relaxation contribution is not small, as it 
changes both the sign and the magnitude of the core-level 
shift on the Mg atom. This happens because the screen- 
ing of the Mg Is core hole is substantially less effective in 
the alloy than in the pure Mg metal. In contrast, the re- 
laxation energy is found to be almost identical in the pure 
metal and the alloy for Au. By inspecting the screen- 
ing density and comparing to free atom calculations, we 
try to offer a simple explanation for the changes in the 
screening properties upon alloying. 

After obtaining a realistic picture of both the initial- 
state and the final-state screening terms, we have tried to 
interpret these results in terms of the potential models. 
However, is has not been possible to reproduce the fea- 
tures of the model. In brief, the potential model assumes 
the following connection between the shift AV in the 
core-level binding energy and the charge transfer Aq: 

AV = (k — M)Aq (1) 

where kAq is the on-site Coulombic potential and — MAq 
is a Madelung term from charges on the other lattice 
sites. Sometimes additional terms are included, e.g. in 
the fonrMta 

= kAq - MAq -6R- e5<j) . (2) 

where —SR is the change in the final-state relaxation en- 
ergy and —S<p is the change in the Fermi energy. Our cal- 
culations can supply reliable values for the well-defined 
quantity SR, but cannot assign values to poorly-defined 
quantities such as 8<j>, k, and M. A substantial effort to 
recover the potential model in either formulation did not 
lead to any quantitative or even qualitative agreement. 



II. CALCULATION AND INTERPRETATION OF 
CORE-LEVEL SHIFTS 

Certain excitation energies (such as core-level shifts 
and atomic ionization energies) can be obtained as the 
difference of the total energies of two self-consistent 
density-functional calculations for the ground state. This 
can be done whenever the excited state is formally the 
ground state for a different set of quantum numbersca 
Although the calculated eigenvalues should not be di- 
rectly associated with the excitation energies, a con- 
nection |-ean be made using Slater's transition state 
concept .E3 

During an experiment such as XPS, an electron is emit- 
ted from the core state into the vacuum. The core-level 
binding energy is the difference of the total energies be- 
tween the unperturbed, homogeneous crystal and the im- 
purity system in which a single atom has a reduced core 
occupation. The first system is easy to handle using stan- 
dard band-structure techniques, whereas the second re- 
quires some treatment suitable for impurities, such as use 
of supercells. 



In a metal, a valence electron moves in from the sur- 
rounding crystal to screen the positive charge of the core 
hole. Effectively, the core electron has been lifted to the 
Fermi level hereby. The energy needed to do this can be 
expected to depend on the position of the core eigenvalue 
before the-excitation and on the degree of screening of the 
core hole.Ej The separation into "initial state" and "final- 
state screening" contributions can be made clearer using 
the transition-state concept. Within DFT, this is done 
using Janak's formula,E3 which states that the deriva- 
tive of the total energy respective to a some occupation 
number equals the corresponding eigenvalue. Applied to 
the present situation, the charge x is taken from the core 
state of one atom in the supercell and put into the valence 
band, and 

=E F - E c {x) = e c (x), (3) 

where Et is the total energy, E-p is the Fermi energy, and 
E c is the core-level eigenvalue. 

Actual calculations show that, to a very good approxi- 
mation, the core eigenvalue drops in a nearly linear fash- 
ion as it is deoccupied, even though the overall core-level 
drop is substantial. For example, e c increases from 1248.6 
to 1364.7 eV when the Mg Is occupation is reduced from 
two to one. Similarly the Au 4/ state starts at 78.6 eV 
below the Fermi energy and drops to 93.3 eV. 

Assuming a strictly linear dependence of e c on x, the 
corel-level binding energy (the change in total energy 
when one electron is taken from the core) can be written 
in various illuminating ways: 

E T (1)-E T (0)= f e c (x)dx (4) 
Jo 

«e c (|) (5) 

«£[e c (0) + e c (l)] (6) 

«e c (0) + |[e c (l)-e c (0)] (7) 

These equations express the full core-level binding en- 
ergy (CLBE) including final-state relaxation effects in 
terms of the eigenvalues at different occupations. Eq. || 
is Slater's transition-state rule, and Eq. || shows that the 
CLBE is the average of the eigenvalues before and after 
removing the core electron. In Eq. |?], the first term e c (0) 
is the initial-state CLBE and the relaxation contribution 
is identified as one-half of the core-eigenvalue drop upon 
depopulation. 

This description is a useful tool to interpret the calcu- 
lated results because it makes contact between core-level 
shifts and differences in the screening response. The core- 
level shift is the difference of the CLBE in two different 
environments, say A and B. The initial-state core- level 
shift is the difference of the static core levels in the un- 
perturbed systems. According to Eq. [5], the final-state 
relaxation contribution can be expressed as the difference 
of the core-eigenvalue drop upon depopulation. In gen- 
eral terms, a core level drops more strongly when the va- 
lence electrons screen the core hole less efficiently. Thus, 
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the initial-state picture for the core-level shift is applica- 
ble if the screening of the core hole is the same in both 
systems. If there is a positive relaxation contribution 
to the core-level shift from A to B, this shows that the 
core-level drop is larger and the screening less effective in 
system B. Conversely, a negative relaxation contribution 
indicates that screening is more effective in B. 



III. CALCULATIONAL PROCEDURE 

To determine the initial-state CLBE, a calculation for 
the unperturbed periodic systems is adequate. Hereby 
it is advantageous to use an all-electron method, which 
gives the core eigenvalues directly. For the complete 
CLBE including final-state relaxation, a supercell is used 
and the difference of the total energy with and without 
a core hole on one atom is evaluated. As discussed, 
the electron taken from the core state is placed into 
the valence band. For additional information, the pro- 
moted charge can take non-integer values. Within a self- 
consistent DFT calculation, the important screening ef- 
fects should be described accurately. Note that the prop- 
erties of the surface, specifically the work function, do 
not enter either description. This must be the case for 
an acceptable model since the true core binding energy 
relative to the Fermi level, expressed as the difference of 
two total energies, is a bulk property. 

The electronic-structure and total-energy calculations 
presented here were done-, with the all-electron full- 
potential LMTO methodjll within the local approxi- 
mation (LDA) to density-functional theoryO Minimiza- 
tion of the energy under a constrained core occupation 
is rigorously justified!!^ in the DFT framework: a self- 
consistent calculation under the chosen constraint pro- 
vides a variational total energy in the parameter subspace 
identified by the constraint. 

We applied this technique to the core levels of the Mg 
Is and Au 4/ levels, first for the pure materials and sec- 
ond for the binary MgAu alloy in the CsCl structure. Ac- 
curate experimental dataffl exist on the core level shifts 
upon formation of this alloy. To make comparisons be- 
tween the materials easier, the fee structure was adopted 
for pure Mg. To study the core-hole-excited solids, we 
used 16-atom supercells for both CsCl-structure MgAu, 
and fee Mg and Au. The distances of the core hole from 
its periodic images exceeds 12 bohr in all cases, and tests 
show that our values for the core level shifts are con- 
verged with respect to cell dimension. The localization 
of the calculated density response to the core-hole per- 
turbation, as detailed below, provides an a posteriori jus- 
tification for the used supercells. The Brillouin-zone in- 
tegration was done using more than 50 irreducible special 
points. Muffin-tin radii for Mg respectively Au are 2.94 
and 2.60 bohr in the pure metals and 2.50 and 2.70 in 
the compound at the experimental lattice constants, and 
are scaled with the lattice constant. All the calculations 



are scalar-relativistic and use the Vosko et al. parame- 
terization of the LDA exchange-correlation potential.LD 



IV. RESULTS FOR MgAu 

The calculated structural parameters for MgAu, Au 
and Mg are given in Table Q, The results for these bulk 
systems are of standard DFT-LDA quality. For the cal- 
culation of the core holes, supercells were built up at the 
theoretical lattice constant. 

Before moving to a discussion of the core level shifts, 
we point out that the absolute core binding energies in 
Table || (referred to Ep, and obtained as total-energy dif- 
ferences) are in remarkable agreement with experiment, 
showing errors below 1%. The data in the Table also il- 
lustrates the large drop of the core eigenvalues when an 
electron is removed (about 14.5 eV for Au and 120 eV for 
Mg). By averaging the eigenvalues before and after re- 
moval of the core electron, Eq. ^ can be verified, showing 
that to a good approximation the core eigenvalue indeed 
drops linearily relative to the Fermi energy as charge is 
removed. 

The calculated initial-state and full core-level shifts are 
obtained by taking the differences of the corresponding 
values in Table ||, leading to the values shown in Table 
III and (graphically) in Figj^. The difference between the 
full and initial-state CLS then gives the screening con- 
tribution. The full results are in good agreement with 
experiment for both cases. We find that the initial-state 
estimate is already accurate for the Au 4/ shift, but that 
it is grossly incorrect and even has the wrong sign for 
Mg Is. The screening contribution to the shift is thus 
completely different for the two types of atom: it is neg- 
ligible for Au, but is the dominant contribution for Mg. 
The screening energies are in reasonable agreement with 
those deduced from Auger parameter measurements, but 
are incompatible with the assumptions made in earlier 
work.Q 

In view of the discussion in Section O, the conclusion 
is that the Au 4/ core hole is screened equally well in 
the pure metal and in the alloy. For Mg, on the other 
hand, the depopulated Is core level has dropped by a 
larger amount in the alloy, showing that the core hole is 
screened significantly less effectively there than in pure 
Mg. Given the size of the effect, an analysis of mea- 
sured core level shifts which does not take screening into 
account is pointless. 

The calculations reproduce the experimental core-level 
shifts and split these unambiguously into an initial-state 
and a final-state screening term. In the rest of this sec- 
tion, we discuss these contributions separately in view of 
the calculated electronic structure and the potential-type 
models. In the end, we will come to the conclusion that 
the potential models are difficult to justify on the basis 
of realistic calculations. 

To help in the interpretation of the results, the site- 
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resolved densities of states are presented in Figs. |2|, |3|, 
and ||. Fig. || compares the electronic structure of Mg, 
Au, and MgAu before making a core hole. The effects of 
a Mg Is or Au 4/ core hole on the valence states of Mg 
and Au are shown in Fig. [j| those in MgAu in Fig. || 



A. Initial-state shifts 

The calculated initial-state core-level binding-energy 
(i-CLBE) shift upon alloying Mg and Au to form MgAu 
is A Mg = -0.45 cV for the Mg Is state and A Au = 0.71 
eV for the Au 4/ state. A charge transfer of about 0. 1-0.2 
electrons from Mg to Au is generally considered reason- 
able, in view of the electronegativity values of 1.31 and 
2.54 for Mg and Au, respectively. A reliable definition 
of the charge transfer (say, as charge density integral) 
from a density-functional calculation is very difficult to 
set up, so we will not try to verify the generally accepted 
value directly. However, we can inspect whether the cal- 
culated initial-state shift is compatible with the potential 
model for this accepted value of the charge transfer. We 
also discuss other modeling concepts which attempt to 
explain the initial-state CLBE shift. In the end, an hon- 
est appraisal is that no simple model can account for the 
calculated values, despite extensive efforts to find one. 

As described above, the basic feature of the poten- 
tial model is that a charge transfer to the atoms of type 
B causes a repulsive on-site Coulomb potential which 
pushes up the core states, reducing the (initial-state) 
CLBE. This effect is only partly compensated by the 
Madelung potential. Thus, the CLBE is reduced for 
those atoms which acquire additional charge, and vice 
versa. On the other hand, we can also present an equally 
simple alternative model, based on a rigid-band descrip- 
tion, where this effect is reversed. Assume (with ref- 
erence to Figs. |[ |[ and [|) that the density of states 
(DOS) of the alloy is obtained by adding together the 
two DOS of the constituents, shifted vertically to line up 
in some way which reflects the bonding. The charge on 
an atom in the alloy is then simply related to the position 
of the shared alloy Fermi energy with respect to the site- 
decomposed valence DOS. Furthermore, we assume that 
the core eigenvalues are at a fixed position relative to the 
valence band, so that the core levels track the shifts of 
the DOS. Since the CLBE is defined relative to the Fermi 
energy, it follows that charge transfer to sites of type B 
is associated with an upward shift of the Fermi energy 
and a larger initial-state CLBE. 

In the case of MgAu, we are assuming that there is 
charge transfer from Mg to Au, and have calculated that 
there is an increase of the Au 4/ CLBE in the alloy. 
The Mg atom has lost some charge by alloying and has 
a reduced CLBE. Even if only the signs are considered, 
these features are incompatible with the potential model, 
but agree with the rigid-band description. In fact, we 
can become ambitious and try use the calculated val- 



ues of the density of states at the Fermi level for the 
pure Mg and Au materials (D_F(Mg)=0.46 states/eV and 
£>f(Au)=0.33 states/eV) to connect the CLBE shifts 
with the charge transfer. Assuming a resonably flat DOS 
at the Fermi energy, the transferred charge for an atom 
of a certain type is approximately equal to 

Aq = D F Ae c 

which yields —0.21 and 0.23 electrons for Mg and Au, 
respectively. These values are close to the generally ac- 
cepted charge transfer for this type of alloy; further- 
more, the charge transferred away from Mg is close to 
the charge transferred to the Au atom. 

Unfortunately, this gratifying result must be consid- 
ered accidental, for several reasons. Foremost is that the 
alloy has a substantially smaller volume than the sum 
of the volumes of the consituents: the cell volumes are 
113.3, 146.9, and 225.9 bohr 3 for Au, Mg, and MgAu, 
respectively. It makes sense to assign the shrinkage of 
13% to the softer Mg atom. Thus, a more correct de- 
scription could be to "prepare" the Mg atom by com- 
pressing it to a smaller volume, then forming the alloy 
from this compressed Mg' and the Au crystal. The total 
i-CLBE shift then is a sum of the effects due to the two 
steps. Independent of whether the shrinkage is assigned 
to the Mg or Au atoms, we are now considering a simpler 
system in which an alloy is formed without any volume 
change. If the rigid-band model is a resonable descrip- 
tion, it should be equally applicable here. When the Mg 
bulk is compressed, the Is core level moves up by 0.51 eV, 
reducing the CLBE from 1248.62 eV to e c =1248.11 eV. 
This value is almost equal to that in the alloy, so that 
the estimated charge transfer for the Mg atom from the 
rigid-band model now comes out close to zero. Unfor- 
tunately, this is not compatible with the charge of 0.23 
electrons added to the Au site, throwing the perceived 
success of the rigid-band model into doubt. 

At this stage, it can be speculated that shifts of the 
Fermi energy should be included, arising from changes 
in the electronic structure due to alloying. Indeed, the 
plots of the DOS in Figs. [| to [| show that the rigid- 
band assumption is not conspicuously well satisfied. To 
demonstrate that all kinds of other effects of similar mag- 
nitude would still be neglected, we focus on just one as- 
pect, namely the role of sp to d charge promotion during 
alloying. This can be most easily investigated in an ASA 
calculation, where the total crystal volume is assigned 
to atomic spheres. The ASA result for the initial-state 
CLBE shifts when alloying Mg' and Au reproduces the 
full-potential calculation reasonably well. The interpre- 
tation of the CLBE shift can now be given a new dimen- 
sion, since we can directly investigate the response of the 
core eigenvalues to changes in the sp and d charges. We 
obtain response parameters de c /dQi « 3 eV/electron for 
the sp and wl.5 eV/electron for the d states. Further- 
more, we can inspect the changes in the partial charges 
Qe when the Mg atom is taken from the pure (com- 
pressed) Mg' crystal and is placed in the alloy, obtaining 
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AQ sp w -0.21 and AQ d w 0.28 electrons. Thus, only 
0.07 electrons are added to the Mg atomic sphere, but 
~ 0.25 electrons are promoted from the sp to the d states. 
Combined with the response parameters, it follows that 
a contribution to the initial-state CLBE shift of about 
—0.3 eV should be attributed to the sp to d promotion. 
Altogether, this shows that not only do the core states 
shift relative to the valence band when charge is trans- 
ferred, but the effect is significantly different depending 
on the angular momentum which takes up the charge. 

In sum, despite attempts in various directions, we have 
not been able to find a simple model which can describe 
the initial-state CLBE shifts caused by alloying. The 
simple potential model is not applicable because even 
the signs are not predicted correctly. The rigid-band 
model seems slightly more plausible, but also suffers 
from a number of shortcomings. Theoretically, an ex- 
tended model could be written down which includes nu- 
merous other relevant effects, such as sp-to-d promotion 
etc. However, this model would be so complicated and 
unwieldy that the overall aim of a simple model would 
be negated. Having confidence in our calculations, we 
believe that the values of —0.45 eV and 0.71 eV for the 
i-CLBE shift are reliable, but we have no convincing way 
to explain these numbers in simple terms. 



B. Final-state screening contribution 

Next, we discuss the final-state screening contribution 
to the CLBE shifts when Mg and Au are alloyed to 
make ordered MgAu. As mentioned above, the calcu- 
lated screening contribution is 0.02 eV for Au and 0.70 
eV for Mg (Table ||). This was interpreted as follows: 
the screening of the Au 4/ core hole happens in a way 
which is nearly independent of the environment. In con- 
trast, the screening of the Mg Is core hole is significantly 
different in the pure Mg bulk and in the alloy. More ex- 
actly, the core hole in the alloy is screened considerably 
less effectively than in pure Mg. In the following, we try 
to analyze this difference in the screening properties. 

A major advantage of an accurate simulation such as 
a DFT calculation is that it can provide data which is 
not accessible to experiment; one example is the separa- 
tion into initial-state and screening contributions. It is 
equally useful to use the calculation as a "microscope" to 
provide quantities such as the DOS or the charge density. 
In the present situation, we can develop a feeling for the 
nature of the core- hole screening by inspecting the screen- 
ing density directly. This is simply the difference of the 
charge density with and without the core hole. One ad- 
ditional electron is in the valence charge, responding to 
the attractive potential of the core hole more or less flex- 
ibly. The screening densities are shown in Figs. ||, [| |^. 
and [v] for pure Mg, pure Au, the Mg core hole in MgAu, 
and the Au core hole in MgAu, respectively. These plots 
are a central result of this paper, making it possible to 



think about the screening cloud in a straightforward and 
unambiguous way. 

By comparing FigsJ^ and ^ for the pure constituents, 
basic differences for Mg and Au are evident. Whereas the 
screening cloud in Mg is wide and extended, screening in 
Au is performed by a localized lump of electrons. For a 
true transition metal, this could be easily explained: the 
screening electron would be taken up by the localized d 
states at the Fermi energy. For Au, however, the d shell 
is already full and this explanation is not possible. 

Instead, the correct explanation for the localised 
screening in Au can be deduced from the corresponding 
DOS plot. On the Au* atom with the core hole, the d 
states are pulled down (and out of the crystal d band) by 
the attractive core-hole potential. The electronic struc- 
ture is similar to that of a Hg impurity in Au. In real 
space, the d states contract, albeit without any change 
in the occupation number. To the screening cloud, this 
process contributes the difference of the contracted and 
uncontracted d shell, which is a positive peak near the 
nucleus surrounded by a negative "ring." At this stage, 
we have not yet taken up the extra screening electron. 
This is done by the sp states which now in turn screen 
(and fill in) the attractive ring. Since the sp states are 
more extended, this cannot be done completely, leaving 
some part of the negative ring visible in the total screen- 
ing cloud. 

For the quality of the screening, the charge closest to 
the nucleus is most relevant. On the Au atom this is 
dominated by the shrinking of the d shell, which can be 
expected to be largely independent of the environment. 
Screening by sp electrons is a more extended affair which 
can be influenced by the environment of the atom. How- 
ever, for Au the sp electrons play a less immediate role, 
even though they actually take up the additional screen- 
ing electron. In contrast, screening of the Mg Is core 
hole is done only by sp valence electrons in the form of 
an extended cloud. Overall, these arguments can explain 
why screening is largely independent of the environment 
in Au but not in Mg, in agreement with the calculated 
results for the alloying process. 

Next, we compare the screening of the core holes in 
the pure materials and the MgAu alloy. For the Au 4/ 
core hole, the screening clouds in the density plots look 
very similar in the central d-electron lump, with some 
differences in the outer regions. Based on the discussion 
above, we can easily accept that the screening is similar 
in Au and MgAu and that only an insignificant contribu- 
tion to the Au 4/ CLBE shift is obtained. For the case 
of Mg, we wish to understand why the screening in the 
alloy is significantly less effective. Unfortunately, in this 
context it is again difficult to obtain a clear answer. 

A first possible explanation for the less effective screen- 
ing in the alloy is that charge has been transferred away 
from the Mg atom, leaving less charge to respond to the 
attractive core hole, leading to reduced screening. How- 
ever, if we count the number of electrons inside a sphere 
of a fixed radius (Rq = 2.8 bohr) we find that the sphere 
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charge in the alloy is 0.36 electrons above that in pure 
Mg. This is presumably a consequence of the reduced 
volume in the alloy. Even though the additional charge 
is mainly in the outer regions of the sphere, it would seem 
to invalidate an explanation based on reduced available 
valence charge. 

Secondly, a comparison of the DOS for pure Mg 
with the Mg site in MgAu shows that the simple-metal 
parabolic sp DOS has changed to some degree of cova- 
lent character in the alloy, with a minimum of the DOS 
around —0.2 Ry. It can be speculated that this leads to a 
somewhat more rigid valence charge density, which can- 
not respond as flexibly to the core hole potential. This 
effect could play a role, but is not confirmed or invali- 
dated by the calculation. 

Finally, Fig.^ shows antiscreening features on the 
neighboring Au atoms. This could be interpreted as 
a "variable wavelength Friedel oscillation" whereby the 
Friedel wavelength changes from a value appropriate to 
Mg to a shorter one on the Au atoms. The antiscreening 
could push the first node of the screening density inwards 
with a corresponding reduction of the screening charge. 
This explanation, while potentially applicable, also can- 
not be confirmed unambiguously. 



V. SUMMARY AND CONCLUSIONS 

In this paper, we have presented the results of ab ini- 
tio density-functional theory calculations of the core level 
shifts which arise upon alloying, using the prototypical 
intermetallic compound MgAu as an example. We were 
interested in the following questions: how well the exper- 
imental results can be reproduced; how the full core level 
shifts can be separated into initial-state and final state 
screening contributions; and whether the results can be 
understood in terms of simple models. 

The agreement to experiment turns out to be good. 
The calculated core-level shifts are 0.73 and 0.25 eV for 
the Au 4/ and Mg Is states, respectively, close to the 
measured values of 0.74 and 0.34 eV. Given the complex- 
ity of the problem, these results are very satisfying. We 
have also found that the absolute core-level binding en- 
ergies, calculated as the difference of two total energies, 
agree to within 1% with the experimental results. 

The calculations give an unambiguous separation of 
each core level shift into a static initial-state contribu- 
tion and a term due to the final-state screening of the 
core hole by the other electrons. Such a separation is cen- 
tral to all subsequent attempts to understand the results 
using simpler concepts. Somewhat unexpectedly, we find 
that the screening contribution is not just a small correc- 
tion, but changes the picture drastically. Specifically, the 
shift of the Mg Is core state changes sign when screening 
effects are included. 

Extensive attempts were made to evaluate the calcu- 
lated results in terms of simpler models. For the initial- 



state shifts, however, no convincing model could be found 
which is able to predict the calculated values. Among 
other considered descriptions, the well-known "simple 
potential model" could not be confirmed. The basic dif- 
ficulty is that a large number of effects influence the core 
level binding energy. The situation is considerably too 
complicated to be cast into any simple model with only 
a few parameters. Possibly, a series of calculations for 
several different systems could uncover trends and help 
to formulate a better model, but we do not consider it 
plausible that an adequate general model can be found. 

For the screening contribution, we have used the calcu- 
lations to obtain accurate images of the screening clouds 
for the different cases. This information cannot be ob- 
tained from experiment and is of major help when trying 
to obtain insight into the nature of the screening process. 
Indeed, straightforward interpretations for the screening 
mechanism at the Mg and Au sites could be deduced. 
Whereas the screening of the Mg core hole is done by 
a relatively extended sp-electron cloud, screening in Au 
takes place in a two-step process. First, the full Au d 
shell contracts in response to the attractive core-hole po- 
tential, then the sp valence electrons fill up the depletion 
ring around the d shell. This description is in line with 
the result that the Mg screening depends on the environ- 
ment, while Au screening does not. 

Two conclusions are drawn from the results. First, 
while it is possible to obtain insight into the electronic 
structure changes upon alloying and the screening be- 
haviour, simple models which try to connect the alloy- 
ing core-level shifts with charge transfer cannot be con- 
firmed. This is mainly due to the complexity of the real 
system which is not compatible with a description involv- 
ing only a few quantities. Specifically, this means that 
charge transfer is only one of several quantities involved, 
and in fact one of the most poorly defined ones. Sec- 
ondly, a full ab-initio calculation can reproduce measured 
core-level binding energies and their shifts to very good 
accuracy. This shows that simpler models are not actu- 
ally needed in order to interpret measured values, where 
such measurements are used to investigate systems with 
unknown properties. Instead, density-functional calcula- 
tions should be used for this purpose. 
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Case 


£c(0) 


«*(!) 


calc. 


exp. 


Au 4/ in Au 


78.64 


93.25 


85.78 


85.88 


Au 4/ in MgAu 


79.35 


94.15 


86.51 


86.62 


Mg Is in Mg 


1248.62 


1364.74 


1306.78 


1303.20 


Mg Is in MgAu 


1248.17 


1365.35 


1307.03 


1303.54 



TABLE II. Calculated and measured core-level binding en- 
ergies in eV, referred to the Fermi energy. Column headings: 
"e c (0)" and "e c (l)" are the calculated eigenvalues before re- 
spectively after making the core hole; "calc" denotes the cal- 
culated CLBE as a total-energy difference; "exp" gives the 
experimental data from Ref. pi 



Case 


initial 


screening 


full 


exp. 


Au 4/ 


0.71 


0.02 


0.73 


0.74 


Mg Is 


-0.45 


0.70 


0.25 


0.34 



TABLE III. Mg Is and Au 4/ core-level shifts in the MgAu 
alloy with respect to pure Au and Mg in eV. Column "initial" 
is the initial-state shift, "screening" is the final-state screen- 
ing contribution, "full" is the full calculation, "exp" is the 
experimental data from Ref. fel 





ao (bohr) 


B (Mbar) 


MgAu th. 


6.09 


1.05 


MgAu exp. 


6.15 




Au th. 


7.68 


1.85 


Au exp. 


7.70 


1.73 


Mg th. 


8.38 


0.40 


Mg exp. 


8.46 





TABLE I. Equilibrium lattice constant and bulk modulus 
for MgAu (CsCl structure), Au (fee) and Mg (fee). The ex- 
perimental Mg lattice constant corresponds to the measured 
volume per atom in the hep structure. 
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Mg with Mg(1s) core hole 



FIG. 1. Graphical representation of the calculated and 
measured core-level binding energy shifts in Table [Ed]. The 
full calculated result in column "Total" is decomposed into 
the intial-state and screening contributions. Note that the 
screening contribution on the Mg atom drastically changes 
the picture. 
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FIG. 3. Calculated site-resolved density of states for Mg 
and Au with, respectively, a Mg Is or Au 4/ core hole. Mg* 
or Au* denotes the site with the core hole, and Inn and 2nn 
are the first and second-nearest-neighbor sites, respectively. 
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FIG. 2. Calculated site-resolved density of states for Mg, 
Au, and MgAu. 
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FIG. 5. Screening charge density for a Is core hole at the 
Mg site in fee Mg. Contour spacing ±5, ±15,. . .times 10 -3 
bohr -3 . The solid and dashed lines show positive and nega- 
tive densities, respectively. 
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FIG. 4. Calculated site-resolved density of states for MgAu 
with a Mg Is or Au 4/ core hole. Details as in Fig.B. 




FIG. 6. Screening charge density for a 4/ core hole at the 
Au site in fee Au. Details as in Fig|B|. 
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